Coulomb Blockade Oscillations in the Thermopower of Open Quantum Dots. 



O 

o 
o 

(N 
X> 

OS 
(N 



A.V. Andreev 1 and K.A. Matveev 2 
1 Department of Physics, University of Colorado, CB 390, Boulder, CO 80309 
2 Department of Physics, Duke University, Durham, NC 27708-0305 
(Draft: February 1, 2008) 

We consider Coulomb blockade oscillations of thermoelectric coefficients of a single electron 
transistor based on a quantum dot strongly coupled to one of the leads. Analytic expression for the 
thermopower as a function of temperature T and the reflection amplitude r in the quantum point 
contact is obtained. Two regimes can be identified: T <C Ec\r\ 2 and T 3> -EcM 2 , where Eq is the 
charging energy of the dot. The former regime is characterized by weak logarithmic dependence of 
the thermopower on the reflection coefficient, in the latter the thermopower is linear in the reflection 
coefficient jr| 2 but depends on temperature only logarithmically. 

PACS numbers: 73.23.Hk, 73.50.Lw 
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Thermoelectric effects in mesoscopic devices have been 
the subject of extensive experimental and theoretical re- 
search The particle- hole asymmetry required for 
such effects can be strongly enhanced in these systems 
as compared to the bulk materials. This and the small 
size of such devices make them promising candidates for 
technological applications, such as micro-refrigerators jj). 
In particular, many experimental and theoretical studies 
in the last few years have focused on the thermoelec- 
tric properties of quantum dots in the Coulomb blockade 
(CB) regime [[|-[| . Most of them concentrated on the CB 
oscillations of the thermopower, S — — where AT is 
the temperature difference across the dot, and AV is the 
voltage necessary to nullify the current. 

The theory of the CB oscillations in the thermopower 
of quantum dots in the weak tunneling regime was con- 
structed in Ref. (J]. This theory takes into account only 
the lowest order tunneling processes, i.e. the sequential 
tunneling, and neglects the cotunneling processes. Its re- 
sults were in agreement with the experiments of Ref. || . 
Later || it became possible to experimentally access 
the regime of lower temperatures and stronger tunnel- 
ing where the cotunneling processes become dominant. 
The theoretical description of this regime was recently 
given in Ref. 0. 

In very interesting recent experiment |8) the CB os- 
cillations in the thermopower of a nearly open quantum 
dot were studied as a function of the reflection coeffi- 
cient |r| 2 in the contact. The setup of these experi- 
ments is schematically represented in Fig. |l|. Surpris- 
ingly, an initial decrease in the amplitude of CB oscilla- 
tions of thermopower with decreasing |r| 2 was followed 
by a plateau with nearly |r|-independent CB oscillations 
of thermopower. This saturation was attributed || to 
the effects of elastic cotunneling || . 

The theory of thermopower for the weak tunneling 
regime developed in Refs. does not apply to this 

case. An additional motivation for studying the thermo- 
electric phenomena in such devices arises from the fact 



that due to the increased transparency of the contacts, 
the open dots are better candidates for micro-refrigerator 
devices Q than the closed ones. 

In this Letter we present a theory of thermoelectric ef- 
fects in a quantum dot in the nearly open regime. We 
consider a quantum dot which is coupled by a tunneling 
junction to the left lead and by a single channel quantum 
point contact (QPC) to the right lead, see Fig. [I]. The 
reflection amplitude in the QPC is assumed to be small, 
|r| <C 1. The mean level spacing 6 in the dot is assumed 
to be vanishingly small. This is a good assumption since 
experimentally |§|] 5 -C T. 

The previous studies of such systems were devoted to 
their thermodynamic and transport properties [^| [L2[ . 
A special feature of the thermoelectric power S is that 
it is sensitive to the average energy transported by elec- 
trons, which in the tunneling approximation depends on 
the odd part of the density of states (DoS) as a function 
of energy. Thus the thermoelectric phenomena represent 
an independent probe of these systems. 

In this Letter we find the thermoelectric coefficient Gt 
of the device in Fig. [j] describing the current response I at 
zero bias, AV = 0, to the difference of the temperatures 
AT between the two leads: Gt — hm |av=o,at^o- 
Our main result is the following expression for Gt of the 
dot: 
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(1) 



Here Gr, <C e 2 /h is the conductance of the left contact, 
e is the absolute value of the electron charge, Eq is the 
charging energy. We have also introduced the energy 
scale T = (8^/n 2 )Ec\r\ 2 cos 2 (nN), which depends on the 
gate voltage N; here In 7 = C w 0.5772 ... is the Euler 
constant. The result was obtained with logarithmic 
accuracy assuming that Ec 3> T, T. 
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FIG. 1. Schematic drawing of a quantum dot electrostati- 
cally defined on a surface of a two-dimensional electron gas. 
The dot is connected to the left lead with temperature T+AT 
by a tunneling junction and to the right lead at temperature 
T by a single channel quantum point contact. The latter is 
characterized by a reflection amplitude r in the contact and 
can be thought of as arising from a point scatterer depicted 
by a cross. The value of r is controlled by the voltage on gate 
R. The voltage on the central gate G determines the optimal 
electron number N in the dot. 

The thermopower S = Gt/G is then obtained from 
Eq. (Q) using the result of Ref. |Q] for the conductance 
G of the device shown in Fig. [l], which we reproduce here 
for completeness: 
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(2) 



For the two limiting cases, T -C T and T ^> V, we 
obtain simplified expressions for the thermopower: 
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for 



t> r, 
T<r. 



(3) 



It is difficult to make a direct comparison of the results 
(0) i (H) with the experiments of Ref. j|] since the experi- 
mental data were presented in terms of a fit to the weak- 
tunneling theory of Ref. Q|. Nevertheless we want to 
point out that even without taking into account the elas- 
tic cotunneling effects || , in the regime T<T the ther- 
mopower only weakly (logarithmically) depends on the 
reflection coefficient, which is consistent with the obser- 
vation in Ref. Q of the thermopower virtually indepen- 
dent of the reflection coefficient. In the opposite regime, 
r <C T, the thermopower is nearly independent of the 
temperature, but scales linearly with the reflection coeffi- 
cient \r\ 2 vanishing, as expected, at perfect transmission. 
Note that even at very low temperatures T <C i^c-lrl 2 , one 
still has T ^> T near half-integer values of the gate volt- 
age N corresponding to the CB peaks of conductance 
The width SN of those regions can be easily found from 
the condition T(N) ~ T. Upon substitution into Eq. 
it gives the estimate of the amplitude of the CB oscilla- 
tions of the thermopower S ~ e ~ 1 \ r \y/T/Ec ln(Ec/T). 

It is interesting to point out that in the low temper- 
ature regime T -C I\ when the conductance (||) shows 



jl2| the temperature dependence G oc T 2 characteris- 
tic of inelastic cotunneling jl3| , the thermopower can be 
expressed in terms of the logarithmic derivative of the 
conductance with respect to the gate voltage 2EcN: 



S = 
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(4) 



This form is analogous to the Cutler-Mott formula |lif| 
for the thermopower of a system of non-interacting elec- 
trons in a metal, but with a different coefficient in front 
of the logarithmic derivative. A similar Cutler-Mott type 
relation holds in the case of weak inelastic cotunneling [Q; 
however, the prefactor of Eq. (|J) contains an additional 
large logarithmic factor ln(Ec/T). In the opposite case 
of high temperature T 3> T no expression similar to the 
Cutler-Mott formula applies. 

Below we present the derivation of the result (Q). 
Following Ref. |l2[ , the electron transport through the 
right QPC can be described by a one-dimensional model 
amenable to bosonization, whereas the left contact can 
be treated in the tunneling approximation. The Hamil- 
tonian of the dot has the form H = Hq + Hr + Hl + He 
where 
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(5a) 
(5b) 

(5c) 
(5d) 



The operators Hr, Hl, and He describe the backscat- 
tcring in the right QPC, tunneling through the left con- 
tact, and the charging energy of the dot, respectively. 
In the equations above a =T,| is the spin label, a pa 
and dka are electron annihilation operators in the dot 
and the left lead respectively, D is the energy cutoff in 
the bosonization, and <p a is the bosonization displace- 
ment operator describing the electron transport through 
the right QPC with II Q being its conjugate momen- 
tum, [<p a (x),Ha'(x')] = iS(x — x')S a , a ' (we have put 
h = 1). The modified form of the tunneling Hamilto- 
nian in Eq. (^) reflects the fact that the electron tun- 
neling event changes the electron number n in the dot. 
This is achieved through the introduction of the charge- 
lowering operator F [0] which satisfies the commutation 
relation [F, n] = F. The current operator through the 
left contact can be obtained from the equation of mo- 
tion for the charge operator / = —era = ie[n,H]. Only 
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Hl contributes to this commutator and gives rise to the 
following expression for the current operator 
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We treat the problem in the lowest order in the tun- 
neling Hamiltonian Eq. (5b). We also assume that the 
conductance of the tunneling contact is much less than 
the conductance quantum, Gl <C e 2 jh. In this approx- 
imation all of the temperature drop happens across the 
left contact. We take the temperature of the left lead to 
be T + AT and that of the dot and the right reservoir to 
be T. In the linear approximation in AT the current I 
can be expressed through the tunneling DoS v(e) in the 
dot as, 
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Here vq is the DoS in the dot in the absence of interac- 



tion, Eq. fl5dp. 

Thus, technically the problem in the tunneling ap- 
proximation reduces to the calculation of the energy- 
dependent tunneling DoS, z/(e). We note that Gt de- 
pends only on the odd (as a function of energy) compo- 
nent of DoS, whereas the conductance G depends only on 
the even one. Therefore, as was mentioned earlier, ther- 
mopower measurements represent an independent test of 
the theory of Coulomb blockade in nearly open dots de- 
veloped in Refs. Moreover, in the leading order 
in max{T, T} / Ec the odd component of the tunneling 
DoS vanishes (TJJ. The thermoelectric coefficient Gt is 
small in the the ratio of max{T, T} /Ec in comparison to 
the conductance G. Its calculation requires going beyond 
the previously adopted approximations |l2j and retaining 
sub-leading order in e/Ec in the tunneling DoS, v(e). 

The tunneling DoS in the dot can be expressed as 
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where Q + itj is the Matsubara Green function, 

Sir) = -J2(Tra pa (r)F(r)al, a (0)F\0)), (9) 



analytically continued to complex time r = j + it. The 



in Eq. 



denote the thermal av- 



angular brackets ( 
erage. 

Because the dynamics of the operators a pa and F 
are decoupled, the Green function in Eq. (|^) factorizes 
into g(r) = G (t)K(t), with G (r) = uqttT/ sin(7rTr) 
being the free electron Green function and K(t) = 
(T T F(r)Ft(0)), @. 



Since the operator F*(0) in K{t) changes the value of 
n from zero to one at t = 0, and F(t) changes it back to 
zero at t = r, the correlator K(t) can be rewritten as 



K(t) = 



Z{r) 

z(Q)' 



(10) 



where Z(t) is a functional integral over 4>uS in the pres- 
ence of the time-dependent charge n T [t) — 6(t)6(T — t). 
Introducing the charge and spin mode variables in the 
right contact (j> C:S (x) — [<t>^(x) ±(j)^(x)\/\/2, we can write 
Z(t) as 



Z(t) 



D[<j> c , cf) s ] exp[-5 c (r) - S , c - 5 , s - S R }. 

(11) 



Here iSo,c+<5o,s represents the free electron part of the ac- 
tion in the absence of backscattering in the QPC, Sc de- 
notes its charging part, and Sr represents the backscat- 
tering in the QPC. These terms are given by 



S c {t) = / dtE, 
Jo 



, 2£>, 
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(12a) 
(12b) 

\/2^ c (0,f)]cos[\/2^ s (0,f)]. (12c) 
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At frequencies below Ec the fluctuations of the charge 
mode, (p c (0,t) are suppressed by the charging energy 



term (12b) and can be integrated out. Furthermore, we 



can evaluate the functional integral over <p c by the saddle 
point ap proximation ignoring the backscattering term, 
Eq. (|l2c| ). 

The action <S sp (r) and the value of the charge mode 
<^ p (0,t) at the saddle point are found to be 
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2jEc sin(7rTY) 



cx '> ' "0.°-"* ^T 
V2(f> s c p (0, t)=ic[N- n T (t)} + T{t) + Tir - t), 
sin(2TmTt) 

„-i n +^T 



m = E 



(13a) 
(13b) 

(13c) 



In Eq. (13a) we have assumed that r ^> E c l , which is 



good approximation since we only need r = 8/2 + it in 
Eq. (§) 



Averaging the backscattering term (12c) over the fluc- 
tuations of d> c we obtain 



8jE c D 
Sr, t = \I- l ^ — \r\ 



dt cos[ V2<f> s c p (0, t)} cos[sf2(j) s (0,t)]. 



Since the charge modes can only be intergated out at fre- 
quencies below the charging energy, one has to assume 
that the energy cutoff in the above action is D ~ Ec- 
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Equation ( |ll| ) can now be written as 

Z{r)=Ne- s ^^Z s {r), 
z s{t) = I D[<j> a ] exp f-<So, s - <Sr )T 



(14a) 
(14b) 



where Af is the r-independent factor which arises from 
the integration over the fluctuations about the saddle 
point and drops out of K(t) in Eq. (jit]). The cor- 
relator K(t) in Eq. jl0| ) then factorizes into K(t) = 
Kq(t)K f (t), where 



- p -5 ap (r) _ 



n 2 T 



2jEc sin(7rTT) 



(15) 



and Kp(f) is the spin part of the correlator which can 
be expressed as 



K f (t) = Z s (t)/Z s {0). 



(16) 



The effective action in Eq. (14b) can be re-fermionized 
following Refs. Jll|Jl^ ]. The Hamiltonian in this repre- 
sentation has the form 

H = iv F [ ^(x)V4>(x)dx + A(t)??[V>(0) - ^(0)], (17a) 



A(t) = -^jv F Ec\r\cos[V2cjfP(t)], 

7T 



(17b) 



where r\ = (c + c') is a Majorana fermion. 

In the limi t T/Ec — * the functions T{t) in 
Eqs. ( 13b| , 13c ) tend to zero, and to the leading or- 
der in T/Ec can be neglected [ Oh Then the time- 
dependent coefficient X(t) in Eq. (|17|) becomes Ao(i) = 
%y/'YVFEc\r\(—l) nT ® cos[ttN]. In this approxima- 
tion |Q the odd component of the tunneling DoS in the 
dot vanishes, thus nullifying the thermopower. There- 
fore we expand Kp(r) in Eq. ( |l6| ) to first order in 
SX(t) = A(t) — Xo(t). In the fermion representation Jl7| ) 
we obtain for the linear in <5A(£) correction to Kp(r) 

f P 

AK f (t)= (-l) n - (t) 5X(t)<S>(T,t)dt, (18a) 
J o 

$(r, t) = (T tV (T)ri(0)r,(t) [^(0, t) - (0, t)]> , (18b) 



where (. . .) denotes the thermal average with the Hamil- 
tonian ( jl7| ) with A = ^y/yvpTi!c\r\cos[TrN] independent 
of time t. 



The average in Eq. ( 18b ) can be evaluated with the 
aid of Wick theorem. It is not difficult to show that 
the thermopower is an odd function of TV. We therefore 
need only to retain the odd in N component A oc jd^_F (t) 
of Eq. (18a). Evaluating the integral in Eq. ( |l8a| ) with 
logarithmic accuracy in Ec /max{T, T} we find: 



A odd K F (r) = -A^^| r | sin(7rjV )ln^ 



£2 + T 2 e K + 1 ■ 



(19) 



The upper energy scale Ec in the logarithmic factor orig- 
inates from the above mentioned energy cutoff D ~ Eq 
of the spin excitations. Using Eq. ( |l9| ) we obtain our 
main result, Eq. (§). 

In conclusion, we have presented a theory of the 
Coulomb blockade oscillations of the thermoelectric co- 
efficient Gt and the thermopower S of quantum dots in 
the anisotropic nearly open regime in the limit where the 
single particle mean level spacing is negligible. Two dis- 
tinct regimes can be identified: the one with T ^> T, and 
the one with r <C T. In the former the thermopower is 
linear in temperature but is nearly independent of the 
reflection coefficient in the QPC and can be expressed 
in the form of Eq. (|J) analogous to the Cutler-Mott for- 
mula E3|. In the latter, the thermopower is linear in the 
reflection coefficient \r\ 2 but depends on the temperature 
only logarithmically. 

We are grateful to B.L. Altshuler, CM. Marcus, 
J.M. Martinis, L.W. Molenkamp and B.Z. Spivak for 
valuable discussions. It is our pleasure to acknowledge 
the warm hospitality of the Aspen Center for Physics, 
the ICTP, Trieste, and the Centre for Advanced Studies 
in Oslo where part of this work was performed. The au- 
thors are A. P. Sloan Research Fellows. A. A. is a Packard 
Research Fellow. This research was supported by the 
NSF Grants No. DMR-9984002 and DMR-9974435. 



[1 

[2 
[3; 
[4; 
[5. 

[6 

[7] 
[8] 

[9 
[10 

[11 
[12 

[13 

[11 



A.V. Anisovich, B.L. Altshuler, A.G. Aronov, and 

A. Yu. Zyuzin, JETP Lett 45, 295 (1987). 

B. Z. Spivak and A.Yu. Zyuzin, Europhys.Lett. 8, 669 
(1989). 

M. Nahum, T.M. Eiles, and J.M. Martinis, Appl. Phys. 
Lett. 65 3123 (1994). 

C. W.J. Beenakker and A. A.M. Staring, Phys. Rev. B46, 
9667 (1992). 

A. A.M. Staring, L.W. Molenkamp, B.W. Alphenhaar, H. 
van Houten, O.J. A. Buyk, M.A.A. Mabesoone, C.W.J. 
Beenakker, and C.T. Foxon, Europhys. Lett. 22, 57 
(1993). 

A.S. Dzurak, C.G. Smith, C.H.W. Barnes, M. Pepper, L. 

Martin-Moreno, C.T. Liang, D.A. Ritchie, and G.A.C. 

Jones, Phys. Rev. B55, R10197 (1997). 

M. Turek and K.A. Matveev, in preparation. 

S. Moller, H. Buhmann, S.F. Godijn, and L.W. 

Molenkamp, Phys. Rev. Lett. 81, 5197 (1998). 

I.L. Aleiner and L.I. Glazman, Phys. Rev. B, 57, 9608 

(1998). 

K. Flensberg, Phys. Rev. B48, 11156 (1993); Physica 

(Amsterdam) 203B, 432 (1994). 

K.A. Matveev, Phys. Rev. B51, 1743 (1995). 

A. Furusaki and K.A. Matveev, Phys. Rev. Lett. 75, 709 

(1995); Phys. Rev. B52, 16676 (1995). 

D. V. Averin and Yu.V. Nazarov, Phys. Rev. Lett. 65, 
2446 (1990). 

M. Cutler and N.F. Mott, Phys. Rev. 181, 1336 (1969). 



